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The Antarctic Muon and Neutrino Detector Array (AMANDA) began collecting data with ten 
strings in 1997. Results from the first year of operation are presented. Neutrinos coming through 
the Earth from the Northern Hemisphere are identified by secondary muons moving upward through 
the array. Cosmic rays in the atmosphere generate a background of downward moving muons, which 
are about 10® times more abundant than the upward moving muons. Over 130 days of exposure, we 
observed a total of about 300 neutrino events. In the same period, a background of 1.05 • 10^ cosmic 
ray muon events was recorded. The observed neutrino flux is consistent with atmospheric neutrino 
predictions. Monte Carlo simulations indicate that 90% of these events lie in the energy range 66 
GeV to 3.4 TeV. The observation of atmospheric neutrinos consistent with expectations establishes 
AMANDA-BIO as a working neutrino telescope. 
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I. INTRODUCTION 
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Energetic cosmic ray particles entering the Earth's at- 
mosphere generate a steady flux of secondary particles 
such as electrons, muons and neutrinos. The electronic 
component of cosmic rays is quickly absorbed. High en- 
ergy muons penetrate the Earth's surface for several kilo- 
meters, while atmospheric neutrinos can easily pass the 
Earth up to very high energies. Interactions of hadronic 
particles, similar to the ones that create the atmospheric 
neutrino flux, will generate neutrinos at sites where cos- 
mic rays are generated and where they interact as they 
travel through the Universe. The goal of observing neu- 
trinos of astrophysical origin determines the design and 
the size of neutrino telescopes. 

The primary channel through which neutrino tele- 
scopes detect neutrinos above energies of a few tens of 
GeV is by observing the Cherenkov light from secondary 
muons produced in f^-nucleon interactions in or near the 
telescope. To ensure that the observed muons are pro- 
duced by neutrinos, the Earth is used as a filter and only 
upward moving muons are selected. A neutrino telescope 
consists of an array of photosensors embedded deeply in a 
transparent medium. The tracks of high energy muons — 
which can travel many hundreds of meters, or even kilo- 
meters, through water or ice — can be reconstructed with 
reasonable precision even with a coarsely instrumented 
detector, provided the medium is sufficiently transparent. 
A location deep below the surface serves to minimize the 
flux of cosmic-ray muons. 

In this paper we demonstrate the observation of at- 
mospheric muon neutrinos with the Antarctic Muon and 
Neutrino Detector Array (AMANDA). These neutrinos 
constitute a convenient flux of fairly well known strength, 
angular distribution, and energy spectrum, which can be 
used to verify the response of the detector. The paper 
will focus on the methods of data analysis and the com- 
parison of observed data with simulations. After a brief 
description of the detector, the data and the methods of 
simulation are introduced in Section III and the general 
methods of event reconstruction are described in Sec- 
IV, Two AMANDA working groups analyzed the 
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data in parallel. The methods and results of both analy- 
ses are described in Sections ^ and VI. Af ter a discussion 
of systematic uncertainties in Section VII we present the 
final results and conclusions. 
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FIG. 1: The present AMANDA detector. This paper de- 
scribes data taken with the ten inner strings shown in ex- 
panded view in the bottom center. 



II. THE AMANDA DETECTOR 

The AMANDA detector uses the 2.8 km thick ice 
sheet at the South Pole as a neutrino target, Cherenkov 
medium and cosmic ray flux attenuator. The detector 
consists of vertical strings of optical modules (OMs) — 
photomultiplier tubes sealed in glass pressure vessels — 
frozen into the ice at depths of 1500-2000 m below the 
surface. Figure ^ shows the current configuration of the 
AMANDA detector. The shallow array, AMANDA-A, 
was deployed at depths of 800 to 1000 m in 1993-94 in 
an exploratory phase of the project. Studies of the op- 
tical properties of the ice carried out with AMANDA- 
A showed a high concentration of air bubbles at these 
depths, leading to strong scattering of light and making 
accurate track reconstruction impossible. Therefore, a 
deeper array of ten strings with 302 OMs was deployed in 
the austral summers of 1995-96 and 1996-97 at depths of 
1500-2000 m. This detector is referred to as AMANDA- 
BIO, and is shown in the center of Fig. |l]. The detector 
was augmented by three additional strings in 1997-98 
and six in 1999-2000, forming the AMANDA-II array. 

In AMANDA BIO, an optical module consists of a sin- 
gle 8" Hamamatsu R5912-2 photomultiplier tube (PMT) 
housed in a glass pressure vessel. The PMT is optically 
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FIG. 2: Variation of tlie optical properties with depth. The 
effective scattering coefficient at a wavelength of 532 nm is 
shown as a function of depth. The z axis is pointing upwards 
and denotes the vertical distance from the origin of the de- 
tector coordinate system located at a depth of 1730 m. The 
shaded areas on the side indicate layers of constant scattering 
coefficient as used in the Monte Carlo. 
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coupled to the glass housing by a transparent gel. Each 
module is connected to electronics on the surface by a 
dedicated electrical cable, which supplies high voltage 
and carries the anode signal of the PMT. For each event, 
the optical module is read out by a peak-sensing ADC 
and a TDC capable of registering up to eight separate 
pulses. The overall precision of measurement of photon 
arrival times is approximately 5 ns. Details of deploy- 
ment, electronics and data acquisition, calibration, and 
the measurements of geometry, timing resolution, and 
the optical properties of the ice can be found in Q . 

The optical properties of the polar ice in which 
AMANDA is embedded have been studied in detail, using 
both light emitters located on the strings and the down- 
going muon flux itself. These studies |^ have shown that 
the ice is not perfectly homogeneous, but rather that it 
can be divided into several horizontal layers which were 
laid down by varying climatological conditions in the 
past Different concentrations of dust in these lay- 
ers lead to a modulation of the scattering and absorption 
lengths of light in the ice, as shown in Fig. |^. The av- 
erage absorption length is about 110 m at a wavelength 
of 400 nm at the depth of the AMANDA-BIO array, and 
the average effective scattering length is approximately 
20 m. 
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FIG. 3: The zenith angle distribution of simulated AMANDA 
triggers per 130.1 days of livetime. The solid line repre- 
sents triggers from downgoing cosmic ray muons generated 
by CORSIKA. The dashed line shows triggers produced by at- 
mospheric neutrinos. 



III. DATA AND SIMULATION 

The data analyzed in this paper were recorded dur- 
ing the austral winter of 1997, from April to November. 
Subtracting downtime for detector maintenance, remov- 
ing runs in which the detector behaved abnormally and 
correcting for deadtime in the data acquisition system, 
the effective livetime was 130.1 days. 

Triggering was done via a majority logic system, which 
demanded that 16 or more OMs report signals within a 
sliding window of 2 /is. When this condition was met, a 
trigger veto was imposed and the entire array read out. 
The raw trigger rate of the array was on average 75 Hz, 
producing a total data set of 1.05 • 10^ events. 

Random noise was observed at a rate of 300 Hz for 
OMs on the inner four strings and 1.5 kHz for tubes on 
the outer six, the difference being due to different levels 
of concentration of radioactive potassium in the pressure 
vessels (details on noise rates can be found in ref Q). A 
typical event has a duration of 4.5 /iS, including the muon 
transit time and the light diffusion times, so random noise 
contributed on average one PMT signal per event. 

Almost all of the events recorded were produced by 
downgoing muons originating in cosmic ray showers. 
Triggers from atmospheric neutrinos contribute only a 
few tens of events per day, a rate small compared to the 
event rate from cosmic ray muons, as shown in Fig. |^. 
The main task of AMANDA data analysis is to seperate 
these neutrino events from the background of cosmic-ray 
muons. Monte Carlo (MC) simulations of the detector re- 
sponse to muons produced by neutrinos or by cosmic rays 
were undertaken to develop techniques of background re- 
jection. 
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Downgoing muons were generated by atmospheric 
shower simulations of isotropic protons with basiev |^ 
or protons and heavier nuclei with CORSIKA using the 
QGSJET generator and tracked to the detector with 
the muon propagation code mudedx Two other 

muon propagation codes were used to check for system- 
atic differences: PROPMU with a 30% lower rate and 
MMC fill with a slightly higher rate. A total of 0.9 • 10^ 
events were simulated. Most characteristics of the events 
generated with basiev were found to be similar to the 
more accurate CDRSIKA-based simulation. For the lat- 
ter, the primary cosmic ray flux as described by Wiebel- 
Sooth and Biermann was used. The curvature of the 
Earth has been implemented in CORSIKA to correctly de- 
scribe the muon flux at large zenith angles. The event 
rate based on this Monte Carlo was 75 Hz and compares 
reasonably well with the observed rate of 100 Hz (after 
deadtime correction). The detector response to muons 
was modeled by calculating the photon fields produced 
by continuous and stochastic muonic energy losses Q, 
and simulating the response of the hardware to these pho- 
tons Upgoing muons were generated by a prop- 
agation of atmospheric neutrinos, which were tracked 
through the Earth and allowed to interact in the ice in 
or around the detector or in the bedrock below p^ . 
Muons that were generated in the bedrock were propa- 
gated using PROPMU until they reached the rock-ice 
boundary at the depth of 2800 m. The muons were then 
propagated through the ice in the same way as those from 
cosmic ray showers. The atmospheric neutrino flux was 
taken from Lipari p^ . 

The Cherenkov photon propagation through the ice 
was modeled to create multidimensional tables of den- 
sity and arrival time probability distributions of the pho- 
ton flux. These photon fields were calculated for pure 
muon tracks and for cascades of charged particles. A 
real muon track was modeled as a superposition of the 
photon fields of a pure muon track and the stochastic 
energy losses based on cascades. The photon fields were 
calculated out to 400 m from the emission point, taking 
into account the orientation of the OM with respect to 
the muon or cascade. In the detector simulation, the ice 
was modeled as 16 discrete layers, as indicated by the 
shaded areas in Fig. ^. The spectral properties of the 
photomultiplier sensitivity, the glass, the gel, and, most 
importantly, the ice itself were included in the simulation 
of the photon propagation. The probability of photon 
detection depends on the Fresnel reflectance at all inter- 
faces, transmittances of various parts, and quantum and 
collection efficiencies of the PMT. The relevant physical 
parameters have been measured in the laboratory, so that 
the spectral sensitivity of the OM could be evaluated. 
Two types of OMs, differing in the type of pressure ves- 
sel, were used in the construction of AMANDA-BIO. The 
inner four strings (AMANDA-B4) use Billings housings 
while the outer six strings use Benthos housings. (Ben- 
thos Inc. and Billings Industries are the manufacturers 
of the glass pressure vessels. Benthos and Billings are 



registered trademarks of the respective companies.) The 
two types of housing have different optical properties. 
The Benthos OMs have an effective quantum efficiency 
of 21% at a wavelength of 395 nm for plane- wave photons 
incident normal to the PMT photocathode. Ninety per- 
cent of the detected photons are in the spectral range of 
345-560 nm. 

An additional sensitivity effect arises from the ice sur- 
rounding the OMs. The deployment of OMs requires 
melting and refreezing of columns of ice, called "hole ice" 
hereafter. This cycle results in the formation of bubbles 
in the vicinity of the modules, which increase scatter- 
ing and affects the sensitivity of the optical modules in 
ways that are not understood in detail. Since the total 
volume of hole ice is small compared to bulk ice in the 
detector (columns of 60 cm diameter, compared to 30 m 
spacing between strings), its effect on optical properties 
can be treated as a correction to the OM angular sen- 
sitivity. The increased scattering of photons in the hole 
ice has been simulated and compared to data taken with 
laser measurements in situ to assess the magnitude of this 
effect. This comparison provides an OM sensitivity cor- 
rection that reduces the relative efficiency in the forward 
direction, but enhances it in the sideways and backward 
directions. The sensitivity in the backward hemisphere 
(90° — 180°) relative to the sensitivity integrated over 
all angles (0° — 180°) of the optical sensor increases from 
20% to 27%, due to this correction, while the average rel- 
ative sensitivity in the forward direction (0° — 90°) drops 
from 80% to 73%. In other words, an OM becomes a 
somewhat more isotropic sensor. 

The effective angular sensitivity of the OMs was also 
assessed using the flux of downgoing atmospheric muons 
as a test beam illuminating both the 295 downward facing 
OMs and the 7 upward facing OMs. We assumed that the 
response of the upward facing OMs to light from down- 
ward muons is equivalent to the response of the down- 
ward facing OMs to light from upward moving muons. 
Based on this assumption we derived a modifled angular 
response function (later referred to as angsens), which 
resulted in a effective reduction of the absolute OM sen- 
sitivity in forward direction. In this model the effective 
relative sensitivity is 67% in the forward hemisphere, and 
33% in the backward hemisphere. This correction will be 
used to estimate the effect of systematic uncertainties in 
the angular response on the final neutrino analysis. 

The simulation of the hardware response included the 
modeling of gains and thresholds and random noise at 
the levels measured for each OM. The transit times of 
the cables and the shapes of the photomultiplier pulses, 
ranging from 170 to 360 ns FWHM, were included in the 
trigger simulation. Multi-photon pulses were simulated 
as superimposed single photoelectron waveforms. In all, 
some 8 • 10^ seconds of cosmic rays were simulated, corre- 
sponding to 7% of the events contained in the 1997 data 
set. 
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IV. EVENT RECONSTRUCTION 



B. Maximum Likelihood Reconstruction 



The reconstruction of muon events in AMANDA is 
done ofHine, in several stages. First, the data are 
"cleaned" by removing unstable PMTs and spurious 
PMT signals (or "hits" ) due to electronic or PMT noise. 
The cleaned events are then passed through a fast filter- 
ing algorithm, which reduces the background of down- 
going muons by one order of magnitude. This reduction 
allows the application of more sophisticated reconstruc- 
tion algorithms to the remaining data set. 

Because of the complexity of the task, and in order to 
increase the robustness of the results, two separate analy- 
ses of the 1997 data set were undertaken. Both proceeded 
along the general lines described above, but differ in the 
details of implementation. The preliminary stages, which 
are very similar in both analyses, are described here. The 
particulars of each analysis will be described in Sections 
0and0. A more detailed description of the reconstruc- 
tion procedure will be published elsewhere ||20|. 



A. Cleaning and Filtering 

The first step in reconstructing events is to clean and 
calibrate the data recorded by the detector. Unstable 
channels (OMs) are identified and removed on a run-to- 
run basis. On average, 260 of the 302 OMs deployed 
are used in the analyses. The recorded times of the hits 
are corrected for delays in the cables leading from the 
OMs to the surface electronics and for the amplitude- 
dependent time required for a pulse to cross the discrim- 
inator threshold. Hits are removed from the event if they 
are identified as being due to instrumental noise, either 
by their low amplitudes or short pulse lengths, or be- 
cause they are isolated in space by more than 80 m and 
time by more than 500 ns from the other hits recorded in 
the event. Pulses with short duration, measured as the 
time over threshold (TOT) , are often related to electronic 
cross-talk in the signal cables or the surface electronics. 
In Analysis II, TOT cuts are applied to individual chan- 
nels beyond the standard cleaning common to both anal- 
yses (see Section VI). 

Following the cleaning and the calibration, a "line fit" 
is calculated for each event. This fit is a simple mini- 
mization of the apparent photon flux direction, for which 
an analytic solution can be calculated quickly |21| (see 
also [||). It contains no details of Cherenkov radiation 
or propagation of light in the ice. Hits arriving at time 
ti at PMT i located at fl are projected onto a line. The 
minimization of = — ~ % ■ U)'^ gives a solu- 

tion for ro and a velocity vif. The results of this fit - at 
the first stage the direction ?/if/|i/if|, at later stages the 
absolute value of the velocity - are used to filter the data 
set. Approximately 80-90% of the data, for which the 
line fit solution is steeply downgoing, are rejected at this 
stage. 



After the data have been passed through the fast fil- 
ter, tracks are reconstructed using a maximum likelihood 
method. The observed photon arrival times do not follow 
a simple Gaussian distribution attributable to electronic 
jitter; instead, a tail of delayed photons is observed. The 
photons can be delayed predominantly by scattering in 
the ice that causes them to travel on paths longer than 
the length of the straight line inclined at the Cherenkov 
angle to the track. Also, photons emitted by scattered 
secondary electrons generated along the track will have 
emission angles other than the muon Cherenkov angle. 
These effects generate a distribution of arrival times with 
a long tail of delayed photons. 

We construct a probability distribution function de- 
scribing the expected distribution of arrival times, and 
calculate the likelihood £timo of a given reconstruction 
hypothesis as the product of the probabilities of the ob- 
served arrival times in each hit OM: 



res I 



ori/ 



(1) 



where ires — ^obs — ^chcr is the time residual (the delay 
of the observed hit time relative to that expected for un- 
scattered propagation of Cherenkov photons emitted by 
the muon), and d± and 6*011 are the distance of the OM 
from the track and the orientation of the module with re- 
spect to the track. The probability distribution function 
p includes the effects of scattering and absorption in the 
bulk ice and in the refrozen ice around the modules. The 
functional form of p is based on a solution to a transport 
equation of the photon fiux from a monochromatic point 
source in a scattering medium p^ , . The free param- 
eters of this function are then fit to the expected time 
profiles that are obtained by a simulation of the pho- 
ton propagation from muons in the ice JT^ , p^ . Varying 
the track parameters of the reconstruction hypothesis, 
we find the maximum of the likelihood function, corre- 
sponding to the best track fit for the event. The result of 
the fit is described by five parameters: three {x,y,z) to 
determine a reference point, and two {9, (j)) for the zenith 
and azimuth of the track direction. Figure ^ shows an 
event display of two upgoing muon events together with 
the reconstructed tracks. 



C. Quality Parameters 

The set of apparently upgoing tracks provided by the 
reconstruction procedure exceeds the expected number 
of upgoing tracks from atmospheric neutrino interactions 
by one to three orders of magnitude, depending on the 
details of the reconstruction algorithm (see Sections 
and VI). In order to reject the large number of "fake 
- events generated by a downgoing muon or 



events 

cascade, but seemingly having an upgoing structure 
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2. Number of Direct Hits, Ndii 

The number of direct hits is defined as the number of 
hits with time delays ires smaller than a certain value. We 
use time intervals of [-15 ns, +25 ns] and [-15 ns,+75 ns], 
and denote the corresponding parameters as N^^^^^ and 



N. 



(75) 



respectively. The negative extent of the window 



allows for jitter in PMT rise times and for small errors in 
geometry and calibration, while the positive side includes 
these effects as well as delays due to scattering of the 
photons. Events with many direct hits (i.e., only slightly 
delayed photons) are likely to be well reconstructed. 



3. Track Length, Ldir 

The track length is defined by projecting each of the 
direct hits onto the reconstructed track, and measuring 
the distance between the first and the last hit. A cut 
on this parameter rejects events with a small lever arm 
for the reconstruction. Direct hits with time residuals 
of [-15 ns,-|-75 ns] are used for the measurement of the 
track length. Cuts on the absolute length, as well as 
zenith angle dependent cuts (which take into account the 
cylindrical shape of the detector) have been used. The 
requirement of a minimum track length corresponds to 
imposing a muon energy threshold. For example, a track 
length of 100 m translates into a muon energy threshold 
of about 25 GeV. 



FIG. 4: Event display of an upgoing muon event. The grey 
scale indicates the flow of time, with early hits at the bottom 
and the latest hits at the top of the array. The arrival times 
match the speed of light. The sizes of the circles corresponds 
to the measured amplitudes. 



we impose additional requirements on the reconstructed 
events to obtain a relatively pure neutrino sample. These 
requirements consist of cuts on observables derived from 
the reconstruction and on topological event parameters. 
Below, we describe the most relevant of the parameters 
used. 



1. Reduced Likelihood, L 

In analogy to a reduced , we define a reduced likeli- 
hood 



L = 



ln£t 



hit 



(2) 



where iVhit — 5, the number of recorded hits in the event 
less the five track fit parameters, is the number of degrees 
of freedom. A smaller L corresponds to a higher quality 
of the fit. 



^. Smoothness, S 

The "smoothness" parameter is a check on the self- 
consistency of the fitted track. It measures the constancy 
of light output along the track. Highly variable apparent 
emission of light usually indicates that the track either 
has been completely misreconstructed or that an under- 
laying muonic Cherenkov light was obscured by a very 
bright stochastic light emission, which usually leads to 
poor reconstruction. The smoothness parameter was in- 
spired by the Kolmogorov-Smirnov test of the consistency 
of two distributions; in our case the consistency of the ob- 
served hit pattern with the hypothesis of constant light 
emission by a muon. 

Figure ^ shows two events to illustrate the character- 
istics of the smoothness parameter. One event is a long 
uniform track, which was well reconstructed. The other 
event is a background event which displays a very poor 
smoothness. 

The simplest definition of the smoothness is given by 



S* = max(|S'j I) where Sj 



1 



N 



i In 



Figure ^ illustrates the smoothness parameter for the two 
events displayed in Fig. ^. Here Ij is the distance along 
the track between the points of closest approach of the 
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FIG. 6: Illustration of the smoothness parameter, which com- 
pares the observed distribution of hits to that predicted for 
a muon emitting Cherenkov light. In a simplest formulation, 
shown here, the prediction is given as a straight line. A large 
deviation from a straight line (0.68) is found for the event on 
the right in Fig. |^. The high quality track-like event on the 
left in Fig. H displays a small deviation (0.09). 



FIG. 5: Two muon events: The upgoing muon event shown 
on the left has a smooth distribution of hits along the track. 
The track-like hit topology of this event can be used to dis- 
tinguish it from background events. The event on the right is 
a background event with a poor smoothness value. 

track to the first and the j"^ hit modules, with the hits 
taken in order of their projected position on the track. 
TV is the total number of hits. Tracks with hits clustered 
at the beginning or end of the track have Sj approaching 
-1-1 or —1, leading to 5 = 1. High quality tracks such as 
the event on the left side of Fig. with S close to zero, 
have hits equally spaced along the track. 



5. Sphericity 

Treating the hit modules as point masses, we can form 
a tensor of inertia for each event, describing the spatial 
distribution of the hits. Diagonalizing the tensor of iner- 
tia yields as eigenvalues li the moments of inertia about 
the principal axes of rotation. For a long, cylindrical 
distribution of hit modules, two moments will be much 
larger than the third. We can reject spherical events, 
such as those produced by muon bremsstrahlung, by re- 
quiring that the normalized magnitude of the smallest 
moment, /i/^/i, be small. 



D. Principal Methods of the Analyses 

The two analyses of the data diverge after the filtering 
stage, following different approaches to event reconstruc- 
tion and background rejection. 

Analysis I uses an improved likelihood function based 
on a more detailed description of the photon re- 
sponse followed by a set of stepwise tightened cuts. 
Analysis 11 uses a Bayesian reconstruction in which 
the likelihood is multiplied by a zenith angle dependent 
prior function, resulting in a strong rejection of downgo- 
ing background. 

Rare backgrounds due to unsimulated instrumental ef- 
fects, such as cross-talk between signal channels and un- 
stable voltage supply, were identified in the course of the 
analyses. These effects either produced spurious triggers, 
or, more often, spurious hits that caused the event to 
be misreconstructed. Different but comparably efficient 
techniques were developed to treat these backgrounds. In 
Analysis I the event topology is inspected; if the spatial 
pattern of hit OMs is inconsistent with the reconstructed 
muon trajectory, the event is rejected. Analysis II at- 
tempts to remove the anomalous hits or triggers through 
identification of characteristic correlations in signal am- 
plitudes and times, which considerably reduces the rate 
of these misreconstructions. 

At this stage the data set in each analysis is reduced to 
several thousand events out of the original 1.05 • 10^, but 
the data are still background dominated. The prediction 
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for atmospheric neutrinos is about 500 at this point. 

For the final selection of a nearly pure sample of 
neutrino induced events, cuts on characteristic observ- 
ables are tightened until the remaining background dis- 
appears. The two analyses use different techniques to 
choose their final cuts, but obtain comparable efficiencies. 
Further details of the analyses can be found in references 
^51, M 



V. ANALYSIS I 

In this analysis the data were processed through three 
levels of initial cuts, designed to reduce the number of 
background events to a manageable size for the final cut 
evaluation. After a first filtering based on the line fit 
(level 1), cuts on the zenith angle, the number of direct 
photons, and the likelihood of the fitted track obtained 
by the maximum likelihood reconstruction were applied 
(level 2). 



Removal of Cascade-Like Events and Detector 
Artifacts 



A third filter level used the results of an iterative like- 
lihood reconstruction with varying track initializations, 
a fit based on the hit probabilities (see Eq. ^ and a re- 
construction to the hypothesis of a high energy cascade, 
e.g. due to a bright seconday muon bremsstrahlung in- 
teraction. 

The first two levels of filtering consisted of relatively 
weak cuts on basic parameters like the zenith angle and 
likelihood. They reduced the data set to about 4 • 10^ 
events. At this stage, residual unsimulated instrumen- 
tal features become apparent, e.g., comparatively high 
amplitude cross-talk produced when a downgoing muon 
emits a bright shower in the center of the detector. Such 
events are predominantly reconstructed as moving verti- 
cally upward and can be identified in the distribution of 
the center of gravity (COG) of hits. Its vertical compo- 
nent (zcog) shows unpredicted peaks in the middle and 
the bottom of the detector (see also Fig. |lj (top), demon- 
strating the effect for Analysis II), while the horizontal 
components (xcoG ^^nd ?/cog) show an enhancement of 
hits towards the outer strings. These strings are read out 
via twisted pair cables, as opposed to the coaxial cables 
used on the inner strings. The twisted pair cables were 
found to be more susceptible to cross-talk signals. Note 
that variations in the optical parameters of the ice due 
to past climatological episodes also produce some vertical 
structure. 

We developed additional COG cuts on the topology of 
the events in order to remove these backgrounds. These 
cuts, which depend on the reconstructed zenith angle, 
use the track lengths L^h and the normalized smallest 
eigenvalues of the tensor of inertia 
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FIG. 7: Characteristic distributions of the center of gravity 
(COG) of events. The figures on the left show the distribution 
of the depth ^cog versus the reconstructed zenith angle. The 
figures on the right show the horizontal location of events in 
the xcoG-2/coG plane of events with m < zcoG < 50 m. 
The positions of the strings are marked by stars. Top: Ex- 
perimental data before application of the COG cuts. Middle: 
Experimental data after application of the COG cuts. Bot- 
tom: Expectation from the BG simulation after cuts. 



Figure shows the different components of the center 
of gravity of the hits and the reconstructed zenith an- 
gle before and after application of the COG cuts, and 
the Monte Carlo prediction for fake upward events stem- 
ming from misreconstructed downgoing muons. The cuts 
remove most of the unsimulated background - in partic- 
ular that far from the horizon - and bring experiment 
and simulation into much better agreement. 

In order to verify the signal passing rates, these cuts 
and those from the previous levels were applied to a sub- 
sample of unfiltered (i.e. downgoing) events but with the 
zenith angle dependence of the cuts reversed, thus using 
the abundant cosmic ray muons as stand-ins for upgoing 
muons. 

In all, these three levels of filtering reduced the data 
set by a factor of approximately 10^ (see Table 



B. Multi-Photoelectron Likelihood and Hit 
Likelihood 

Before the final cut optimization the last, most elab- 
orate reconstruction was applied, combining the likeli- 
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hoods for the arrival time of the first of muUple photons 
in a PMT with the hkehhoods for PMTs to have been 
hit or have not been hit. 

T he pro babihty densities p(tres | d'f', 6'^^.-) (see eq.l, Sec- 



tion EVB) describe only the arrival times of single pho- 
tons. Density functions for the multi-photoelectron case 
have to include the effect of repeatedly sampling the dis- 
tribution of photon arrival times. For several detected 
photons, the first of them is usually less scattered than 
the average photon (which defines the single photoelec- 
tron case). Therefore the leading edge of a PMT pulse 
composed of multiple photoelectrons (MPE) will be sys- 
tematically shifted to earlier times compared to a single 
photoelectron. The MPE likelihood C^^fJ^ [|| uses the 
recorded amplitude information to model this shift. 

In the reconstructions mentioned so far, the timing in- 
formation from hit PMTs was used. However, a PMT 
which was not hit also delivers information. The hit like- 
lihood Chit does not depend on the arrival times but rep- 
resents the probability that the track produced the ob- 
served hit pattern. It is constructed from the probability 
densities Phit(d^\ ^^orl) ^^^^ ^ given PMT i was hit if it 
was in fact hit, and the probabilities (1 — Phit(c?x\ 6*01!) 
that a given PMT j was not hit if it was not hit: 



Chit — 



Not. 



>W d'h \ (4) 



where the first product runs over all hit PMTs and the 
second over all non-hit PMTs. 

The likelihood combining these two probabilities is 



r — /'MPE 

•'- ~ •'-time ■ ■'-hit 



(5) 



A cut on the reconstructed zenith angle obtained from 
fitting with C leaves less than 10^ events in the data set, 
defined as level 4 in Fig. |[ 



C. Final Separation of the Neutrino Sample 

For the final stage of filtering, a method (CutEval) 
was developed to select and optimize the cuts taking into 
account correlations between the cut parameters. A de- 
tailed description of this method can be found in [ p7[ . 
The principle of CutEval is to numerically optimize the 
ratio of signal to -^/background by variation of the selec- 
tion of cut parameters, as well as the actual cut values. 
Parameters are used only if they improve the efficiency 
of separation over optimized cuts on all other already in- 
cluded parameters. A first optimization was based purely 
on Monte Carlo, with simulated atmospheric neutrinos 
for signal and simulated downgoing muons forming the 
background. This optimization yielded four such inde- 
pendent parameters. Two other optimizations involved 
experimental data. In both cases, experimental data have 



been defined as the background sample. In one case, the 
signal was represented by atmospheric neutrino Monte 
Carlo, in the other by experimental data subjected to 
zenith angle inverted cuts (i.e. to downward events pass- 
ing the quality cuts, but being "good" events with respect 
to the upper hemisphere instead - like neutrino candi- 
dates - with respect to the lower hemisphere). These 
latter optimizations yielded two additional parameters, 
which rejected a small contribution of residual unsimu- 
lated backgrounds: coincident muons from simultaneous 
independent air showers and events accompanied by in- 
strumental artifacts such as cross-talk. After application 
of these two cuts to simulated and experimental data, 
the distributions of observables agree to a satisfactory 
precision. 

Once the minimal set of parameters is found, the op- 
timal cut values can be represented as a function of the 
number of background events A'bg passing the cuts. The 
result is a path through the cut parameter space which 
yields the best signal efficiency for any desired purity of 
the signal, characterized by Nbg- Using this represen- 
tation, one can calculate the number of events passing 
the cuts as a function of the fitted Nbg for signal and 
for background Monte Carlo. Figure || (top) shows this 
dependence for simulations as well as for experimental 
data, with A^bg varying from trigger level to a level that 
leaves only a few events in the data set. One observes 
that the actual background expectation falls roughly lin- 
early as the fitted A^bg is reduced. Below values of a few 
hundred events the signal is expected to dominate the 
event sample. The experimental curve follows the ex- 
pectation from the sum of background and signal Monte 
Carlo. For large A^bGj the observed event rate follows 
the background expectation. At smaller A'bgj the ex- 
perimental shape turns over into the signal expectation 
and follows it nicely down to the sample of events with 
highest quality (the smallest values of A^bg)- For a mod- 
crate background contamination of Abg = 10; one gets 
a total of 223 neutrino candidates. The parameters and 
cut values as obtained by the CutEval procedure are 
summarized in Table ^. 

Figure || (bottom) translates the background parame- 
ter A^BG into an event quality parameter Q, defined as 
Q = ln{No/NBG) = ln(1.05 • IOVA^bg). The plot 
shows the ratios of events from the upper figure as a 
function of Q. At higher qualities {Q > 17), the ratio of 
observed events to the atmospheric neutrino simulation 
flattens out with a further variation of only 30%. The 
value at Q — 17 is approximately unity for the angsens 
Monte Carlo and about 0.6 for the standard Monte Carlo 
(chosen in Fig. |^, top) and approximately unity for the 
angsens Monte Carlo (chosen in Fig. ||, bottom). 

Table |l| lists the cut efficiencies for the atmospheric 
neutrino simulation (with and without the implementa- 
tion of the angular sensitivity fitted model angsens of 
the OMs — see Sections III and VII), the background 



simulation of atmospheric muons from air showers (with- 
out angsens) and the experimental data. Again, the 
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Parameter 


Cut 


Explanation 


\s\ 


< 0.28 


see Section IV C 4 


1'-' 1/ \^"mpe -^w y 


< 0.01 


Tightens the requirement on the smoothness for tracks 
close to the horizon where background is high. 


(A^dir - 2) ■ Ldir 


> 750m 


Lever arm of the track times the number of supporting points. 


lo? / /"^ 1 


< -7.7 


Ratio of the likelihoods of the best 
upgoing and best downgoing hypotheses. 


*(mpe,lf) 


< 35° 


Space angle between the results from the multi-photon 
likelihood reconstruction and the line fit. This cut 
effectively removes cross-talk features. 




< 0.55 


Parameter combining the two smoothness definitions (here 
calculated using only direct hits). This cut effectively removes 
coincident muon events from independent air showers. 



TABLE I: Final quality parameters and cuts obtained from the cut evaluation procedure. The "direct" time interval for variables 
NdiT, Ldir, and Sdir is [-15 ns,-|-75 ns]. The first four rows show cut parameters obtained by all (Monte Carlo and experimental) 
searches, the last two rows show two additional (weaker) cuts, which were found to remove unsimulated backgrounds. 



Filter level 


Atm. u 


Atm. f MC 


Atm. fi MC 


Experimental 




MC 


angsens 


(Background) 


data 


Events at trigger level 


8978 


5759 


9.03 ■ lO** 


1.05 ■ lO'-" 


Efficiency at Level 1 


0.34 


0.37 


0.4 TO"' 


0.5 TO"' 


Efficiency at Level 2 


0.15 


0.15 


0.4 TO"^ 


0.4 TO"^ 


Efficiency at Level 3 


0.7 TO^^ 


0.7 TO"^ 


0.7 TO^^ 


0.1 TO"* 


Efficiency after final cuts 


0.4 TO'^ 


0.4 40"^ 


0.6 TO-** 


0.2 TO"'' 


No. of events 


362 ±4 


237 ±6 


7±5 


223 


passing final cuts 






normalized 





TABLE II: The cut efficiencies for the atmospheric neutrino Monte Carlo (MC), the atmospheric muon background Monte 
Carlo, and the experimental data for 130 days of detector livetime. Efficiencies are given for filter levels LI to L4. L4 is the 
final selection. All errors are purely statistical. The final background prediction of 7 events has been normalized at trigger 
level. 



experimental numbers agree well with the background 
simulation up to the first two filter levels. Later, the 
Monte Carlo underestimates the experimental passing 
rates slightly. The last row shows the expected num- 
bers of events for the last stage of filtering. If, in addi- 



tion, the effect of neutrino oscillations (see Section VII ) is 
included, the atmospheric neutrino simulation including 
the angsens model predicts 224 events, in closest agree- 
ment with the experiment. However, the 5% effect due 
to oscillations is smaller than our systematic uncertainty 
(see Section VII). 



D. Characteristics of the Neutrino Candidates 



Time Distribution 



downward muon rate during the Antarctic winter Pq| and 
should result in a 10% deficit of triggers with respect to 
upward neutrino events.) 



2. Zenith Angle Distribution 

Figure [ic| shows the zenith angle distribution of the 223 
neutrino candidates compared to the Monte Carlo predic- 
tion for atmospheric neutrinos and the few remain- 
ing events predicted by background simulations. Note 
that Monte Carlo is normalized to experiment. (The to- 
tal number of events is 362 for the atmospheric neutrino 
simulation and 223 for experiment, i.e. there is a deficit 
of 39 percent in the absolute number of events.) There 
is good agreement between the prediction and the exper- 
iment in the shape of the angular distribution. 



Figure ^ shows the cumulative number of neutrino 
events as well as the cumulative number of event trig- 
gers plotted versus the day number in 1997. One can 
observe that the neutrino events follow the number of 
triggers, albeit with a small deficit during the Antarctic 
winter. This deficit is consistent with statistical fiuctua- 
tions. (Actually, seasonal variations slightly decrease the 



3. Characteristic Distributions and Visual Inspection 

Four methods were used to evaluate the effectiveness 
of the analysis and the level of residual backgrounds: 
a) A'^— 1 cuts^ b) unbiased variables, c) low level distribu- 
tions, and d) visual inspection. 
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FIG. 8: The fitted background parameter Nbg- Top; Num- 
ber of events versus TVbg- Smaller values of A^bg correspond 
to harder cuts. Below A^bg = 1500 the CutEval parameter- 
ization was used to calculate the cut values corresponding to 
A^BG • For larger values of A^bg the data points correspond to 
the cuts from the filter levels: Level 4 (see Section V B), Level 
3, Level 2, Level 1, and trigger level (table Bottom: Ra- 
tios of events passing in the experimental data compared to 
various Monte Carlo expectations for signal and background 
as a function of event quality. The dashed line indicates the 
final cuts. 



FIG. 9: The integrated exposure of the AMANDA detector 
in 1997. The figure shows the cumulative number of triggers 
(upper curve) and the number of observed neutrino events 
(lower curve) versus the day number. The intervals with zero 
gradient correspond to periods where the detector was not 
operating stably; data from these periods were excluded from 
the analysis. 




> 



* Lxp 

□ Sicnal fvIC 
O BG MC 




FIG. 10: Zenith angle distribution of the experimental data 
compared to simulated atmospheric neutrinos and a simulated 
background of downgoing muons produced by cosmic rays. 
In this figure the Monte Carlo prediction is normalized to 
the experimental data. The error bars report only statistical 
errors. 



a) The N — 1 test evaluates the N final cuts one 
by one and yields an estimate of the background con- 
tamination in the final sample. One applies all but one 
of the final cuts (the one in the selected variable), and 
plots the data in this variable. In the signal region of 
this variable (defined by the later applied cut) shapes 
of experiment and signal Monte Carlo should agree. In 



the background region, the experimental data should 
approach the expected background shape. Figure |ll| 
shows four of these distributions. The applied cut is 
shown by a dotted line. All four cuts satisfy the test: 
the shape of the distributions agree reasonably well on 
both sides of the applied cuts. Two — 1 distribution 
from Analysis II are shown in Fig. n3. 
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FIG. 11: Two distributions of variables used as cut param- 
eters in the last filter level (see Table 1 for an explanation of 
the variables). In both cases, all final cuts with the excep- 
tion of the variable plotted have been applied. The cuts on 
the displayed parameters are indicated by the deished vertical 
lines. Arrows indicate the accepted parameter space. 
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FIG. 12: Distributions of zcoG for the experiment and atmo- 
spheric neutrino signal Monte Carlo. MC standard and MC 
bulk ice denote two different ice models. The first includes 
vertical ice layers in accordance with Fig. 2, the second uses 
homogeneous ice. 



However, as can bee seen from Fig. gO[ part of these 
disagreements disappear if the the standard atmospheric 
neutrino MC is replaced by the angsens MC version. 

c) In order to account for possible pathological low 
level features in the data sample (especially cross-talk), 
we i) investigated basic pulse amplitude and pulse 
width (TOT) distributions and ii) re- fitted all events 
after the cross-talk hit cleaning procedure applied in 
Analysis II (which is tighter than the s tandard cross-talk 
cleaning introduced in Section IV A ). Both: 



these 



ated zenith angles 



distributions and that for the recalcu 
show no significant deviation from the previous ones. No 
cross-talk features are found in the resulting neutrino 
sample. 

d) Finally, a visual inspection of the full neutrino 
sample was performed, by visually displaying each event 
like in Fig. 4. The visual inspection gives consistent 
results with the other methods of background estimation 
and yields an upper limit on the background contamina- 
tion of muons from random coincident air showers (see 
below) . 



b) An obvious test is the investigation of distribu- 
tions of unbiased variables (i.e. variables to which no 
cuts have been applied) in the final neutrino sample. 
Here, the experimental distributions follow the Monte 
Carlo signal expectations nicely. Some deviations are 
observed, especially in the number of OMs hit and the 
velocity v\{ obtained from the line fit (see Section IV, A). 



E. Background estimation 



The results of four independent methods of background 



estimation are summarized in tabic HI 



First, the background Monte Carlo itself gives an es- 
timate. It yields 7 events if rates are normalized to the 
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trigger level (see table ||). Because the passing rates 
differ slightly between the experiment (higher) and the 
background Monte Carlo (lower) , we made the conserva- 
tive choice to renormalize the background Monte Carlo 
to the level 3 experimental passing rate. This gives an 
estimate of about 16 background events in our final sam- 
ple. 



BG estimation method 


estimation 


BG MC 


16 


±8 


N-1 cuts 


14 


±4 


zcoG distributions 


< 


35 


Visual inspection 


< 


23 



TABLE III: Various estimates of the background remaining 
in the experimental data sample of 223 neutrino candidates. 



VI. ANALYSIS II 

The second analysis follows a different approach; in- 
stead of optimizing cuts to reject misreconstructed cos- 
mic ray muons, this analysis concentrates on improv- 
ing the reconstruction algorithm with respect to back- 
ground rejection. The large downgoing muon flux implies 
that even a small fraction of downgoing muons misrecon- 
structed as upgoing will produce a very large background 
rate. Equivalently, for each apparently upgoing event, 
there were many more downgoing muons passing the de- 
tector than there were upgoing muons; even though any 
single downgoing muon had only a small probability of 
faking an upgoing event, the total probability that the 
event was a fake is quite high. 



A. Bayesian Reconstruction 



From the iV — 1 distributions we obtained an alterna- 
tive approximation of the residual background. We re- 
normalized both signal and background MCs in the back- 
ground region to fit the number of experimental events 
in the background region. The number of re-normalized 
background MC events in the signal region is then a back- 
ground estimate. This estimate was performed N times 
(once for each iV — 1 distribution) . The average over all 
N estimations yields 14 background events. Note that 
this averaging procedure is reasonable only for the case 
of independent cuts. With the method by which we have 
chosen the cut parameters, this condition is fulfilled to 
first approximation. 

We have found that cross-talk hits are related to the 
characteristic triple-peak structure in the distribution of 
the vertical component of the center of grav ity o f hits 
(zcog) which has been discussed in Section VA - see 



Fig. 7 and also Fig.y^ (top). Since there are remaining 
cross talk hi ts whi ch have survived the standard cleaning 
(see section IV A), this distribution was studied in de- 
tail. As shown in Fig. the final experimental sample 
of neutrino candidates shows no statistically significant 
excess with respect to the atmospheric neutrino Monte 
Carlo prediction in the regions of the characteristic peaks. 
Therefore, an upper limit on this special class of back- 
ground was derived and yields < 35 events. 

The visual inspection of the neutrino sample yields 13 
events. Seven of them show the signature of coincident 
muons from independent air showers; i.e., two well sep- 
arated spatial concentrations of hits, each with a down- 
ward time flow but with the lower group appearing earlier 
than the upper one. Taking into account the scanning ef- 
ficiencies which were determined by scanning signal and 
background Monte Carlo events, an upper limit of 23 
events is obtained from visual inspection. 

Combining the results from the above methods, the 
expected background is estimated to amount to 4 to 10% 
of the 223 experimental events. 



This analysis of the problem motivates a Bayesian ap- 
proach to event reconstruction. Bayes' Theorem in 
probability theory states that for two assertions A and 
B, 

P{A I B) P{B) = P{B I A) P{A), 

where P{A \ B) is the probability of assertion A given 
that B is true. Identifying A with a particular muon 
track hypothesis /i and B with the data recorded for an 
event in the detector, we have 

P(^ I data) = /:timc(data | ^) P{p), 

where we have dropped a normalization factor P(data) 
which is a constant for the observed event. The function 
'Ctimc is the regular likelihood function of Eq. ^ and P(/i) 
is the so-called prior function, the probability of a muon 
/I = ^(x, y, z, 9, (j>) passing through the detector. 

For this analysis, we have used a simple one- 
dimensional prior function, containing the zenith angle 
information at trigger level in Fig. ^. By accounting in 
the reconstruction for the fact that the flux of downgo- 
ing muons from cosmic rays is many orders of magnitude 
larger than that of upgoing neutrino-induced muons, the 
number of downgoing muons that are misreconstructed 
as upgoing is greatly reduced. It should be noted that 
the objections that are often raised with respect to the 
use of Bayesian statistics in physics are not relevant to 
this problem: the prior function is well defined and nor- 
malized and independently known to relatively good pre- 
cision, consisting only of the fluxes of cosmic ray muons 
and atmospheric muon neutrinos. 



B. Removal of Instrumental Artifacts 

The Bayesian reconstruction algorithm is highly effi- 
cient at rejecting downgoing muon events. Of 2.6 • 10^ 
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events passing the fast filter, only 5.8 • 10^ are recon- 
structed as upgoing. By contrast, the standard maxi- 
mum likelihood reconstruction produces about 2.4 • 10^ 
false upgoing reconstructions. However, less than a thou- 
sand neutrino events are predicted by Monte Carlo, so it 
is clear that a significant number of misrcconstructions 
remain. 

Detailed inspection of the 5.8 • 10'* events reveals that 
the vast majority is produced by cross-talk overlaid on 
triggers from downgoing muons emitting bright stochas- 
tic light near the detector. This cross-talk confuses the 
reconstruction algorithm, producing apparently upgoing 
tracks. Because cross-talk is not included in the detec- 
tor simulation, the characteristics of the fakes are not 
predicted well by the simulation, and the rate of misre- 
construction is much higher than predicted. 

The cross-talk is removed by additional hit cleaning 
routines developed by examination of this cross-talk en- 
riched data set. For example, cross-talk in many channels 
can be identified in scatter plots of pulse width vs. ampli- 
tude, as shown in Fig. |l^. The pulse width is measured 
as time-over-threshold (TOT). Real hits form the distri- 
bution shown on the left. High amplitude pulses should 
have large pulse width. This is not the case for cross-talk 
induced pulses. In channels with high levels of cross-talk, 
an additional vertical band is found at high amplitudes 
but short pulse widths, as seen in the lower figure. 

Other hit cleaning algorithms use the time correlation 
and amplitude relationship between real and cross-talk 
pulses and a map of channels susceptible to cross-talk and 
the channels to which they are coupled. An additional 
instrumental effect, believed to be caused by fluctuating 
high voltage levels, produces triggers with signals from 
most OMs on the outer strings but none on the inner 
four strings; some 500 of these bogus triggers were also 
removed from the data set. The 5.8 • 10* upgoing events 
were again reconstructed after the additional hit cleaning 
was apphed. Only 4.9- lO'^ (8.4%) of the events remained 
upgoing, compared to an expectation from Monte Carlo 
of 1855 atmospheric muon events (37.8% of the total be- 
fore the additional cleaning), and 555 atmospheric neu- 
trino events. Figure |lj (top) shows that while there has 
been a significant reduction in the instrumental back- 
grounds, an unsimulated structure still remains in the 
center-of-gravity (COG) distribution for these remaining 
data events. The application of additional quality cri- 
teria brings this distribution in agreement, as shown in 
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Figure 14 (bottom). 



C. Quality Cuts 

The improvements in the reliability of the reconstruc- 
tion algorithm described above obviated the need for 
large numbers of cut parameters or for careful opti- 
mization of the cuts. Because the signal-to-noise of the 
upward-reconstructed data is quite high to begin with, 
we have the possibility of comparing the behavior of real 
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FIG. 13: Pulse amplitude vs. duration for modules on the 
outer strings. Normal hits lie in the distribution shown in 
the upper figure. High amplitude pulses of more than a few 
photoelectrons are valid only if the pulse width is also large. 
Cross-talk induced pulses of high amplitude are characterized 
by small time-over-threshold (TOT) . The cut-off seen at high 
amplitude is due to saturation of the amplitude readout elec- 
tronics. 



and simulated data over a wide range of cut strengths 
to verify that the data agree with the predictions for up- 
going neutrino-induced muons, not only in number but 
also in their characteristics. Using the cut parameters 



described in Section IV C (with the likelihood replaced 
by the Bayesian posterior probability) and a requirement 
that events fitted as relatively horizontal by the line fit fil- 
tering algorithm not be reconstructed as steeply upgoing 
by the full reconstruction (a requirement that suppresses 
residual cross-talk misrcconstructions) , an index of event 
quality was formed. 

To do so, we rescale the six quality parameters de- 
scribed above by the cumulative distributions of the sim- 
ulated atmospheric neutrino signal, and consider the six- 
dimensional cut space formed by the rescaled parameters. 
A point in this space corresponds to fixed values of the 
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FIG. 14; Top: Event center of gravity distribution after recon- 
struction with special cross-talk cleaning algorithms applied 
to the events. Unsimulated background remains. Bottom: 
The data agree with the neutrino signal after application of 
additional quality cuts. 



quality parameters, and events can be assigned to loca- 
tions based on their track length, sphericity, and so forth. 

It is difficult to compare the distributions of data and 
simulated up- and downgoing muons directly because 
of the high dimensionality of the space. We therefore 
project the space down to a single - "quality" - dimen- 
sion by dividing it into concentric rectangular shells, as 
illustrated in Fig. The vertex of each shell lies on a 
line from the origin through a reference set of cuts which 
are believed to isolate a fairly pure set of neutrino events. 
Events in the full cut space are assigned an overall quality 
value, based on the shell in which they lie. 

With this formulation we can compare the character- 
istics of the data to simulated neutrino and cosmic-ray 
muon events. Figure ^ compares the number of events 
passing various levels of cuts; i.e., the integral number 
of events above a given quality. At low qualities, g < 3, 
the data set is dominated by misreconstructed downgoing 
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FIG. 15: Definition of event quality. Events are plotted in 
A^- dimensional cut space (two dimensions are shown here for 
clarity). A line is drawn from the origin (no cuts) through 
a selected set of cuts, and the space is divided into rectan- 
gular shells of equal width. Events are assigned a quality q 
according to the shell in which they are found. 



muons, data as well as the simulated background exceed 
the predicted neutrino signal. At higher qualities, the 
passing rates of data closely track the simulated neutrino 
events, and the predicted background contamination is 
very low. 
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FIG. 16: Numbers of events above a certain quality level, for 
downgoing muon Monte Carlo, atmospheric neutrino Monte 
Carlo, and experimental data. 

We can investigate the agreement between data and 
Monte Carlo more systematically by comparing the dif- 
ferential number of events within individual shells, rather 
than the total number of events passing various levels of 
cuts. This is done in Fig. ^ where the ratios of the 
number of events observed to those predicted from the 
combined signal and background simulations are shown. 
One can see that at low quality levels there is an excess 
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FIG. 17: Ratio of data to Monte Carlo (cosmic ray muons plus 
atmospheric neutrinos) . Unlike Fig. 1^, the plot is differential 
— the ratio at a particular quality does not include events at 
higher or lower qualities. 



in the number of misreconstructed events observed. This 
is mainly due to remaining cross-talk. There is also an 
excess, though statistically less significant, at very high 
quality levels, which is believed to be caused by slight 
inaccuracies in the description of the optical parameters 
of the ice. Nevertheless, over the bulk of the range there 
is close agreement between the data and the simulation, 
apart from an overall normalization factor of 0.58. The 
absolute agreement is consistent with the systematic un- 
certainties. It should be emphasized that the quality pa- 
rameter is a convolution of all six quality parameters, and 
so the flat line in Fig. ^ demonstrates agreement in the 
correlations between cut parameters. 



D. Background Estimation and Signal Description 

If we reduce the 4,917 upward-reconstructed events by 
requiring a quality of at least 7 on the scale of Fig. |l^, we 
obtain a set of 204 neutrino candidates. The background 
contamination, which is due to misreconstructed down- 
going muons, was estimated in three ways. The first way 
is to simulate the downgoing muon flux, bearing in mind 
that we are looking at a very low tail (10~*) of the to- 
tal muon distribution. The second way is to renormalize 
the signal simulation by the factor of 0.58 obtained from 
Fig. |l^ and subtract the predicted events from the ob- 
served data set (accepting the excess at extremely high 
qualities, however, as signal). The third way, a cross 
check on the first two methods, is to examine the data 
looking for fakes due to unsimulated effects such as cross- 
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FIG. 18: The zenith angle distribution of upward recon- 
structed events. The size of the hatched boxes indicates the 
statistical precision of the atmospheric neutrino simulation. 
The Monte Carlo prediction is normalized to the data. 



talk, independent coincident downgoing muons, and so 
forth. All three methods yield estimates of 5-10% con- 
tamination. 

The zenith angle distribution for the 204 events is 
shown in Figure ^_8|, and compared to that for the simu- 
lation of atmospheric neutrinos. In the figure the Monte 
Carlo events are normalized to the number of observed 
events to facilitate comparison of the shapes of the distri- 
butions. The agreement in absolute number is consistent 
with the systematic uncertainties in the absolute sensi- 
tivity and the flux of high energy atmospheric neutrinos. 
The shape of the distribution of data is statistically con- 
sistent with the prediction from atmospheric neutrinos. 
Figure |l^ (bottom) shows the distribution of the zcoG 
parameter for the 204 events. The level 7 quality cuts 
have removed the remainder of the instrumental events 
left after the Bayesian reconstruction with the improved 
cross-talk cleaning algorithm, bringing the data events 
in line with the atmospheric neutrino expectations. The 
efficiencies corresponding to the three steps of the data 
analysis: 1) events reconstructed upward, 2) events re- 
constructed upward with cross-talk cleaning, and 3) with 
additional level 7 quality cuts are summarized in table 

ivl. 



Figure 19 (top) shows the smoothness distribution for 
events that have passed the quality level 7 cuts for the 
five observables except smoothness. The vertical dashed 
line at smoothness ^ 0.29 shows the value of the level 7 
smoothness cut. This cut removes the tail of fake events 
leaving a good agreement between remaining data and 
Monte Carlo. Figure |l^ (bottom) shows the same plot 
for the direct length variable. Again, a clear tail of fake 
events is removed by requiring a direct length of greater 
than 70 meters. 
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8978 
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Erfaciency: Reconstructed upgomg 


0.55 • 10"^ 




0.55 ■ 10""' 


Efficiency: Reconstructed upgoing 
(w/ cross-talk cleaning) 


(2.1 ±0.08) ■ lO"** 


(6.2 ±0.06) • 10"^ 


4.7 ■ 10"'^ 


Efficiency: Final Cuts (q>7) 


(1.9 ±0.6) ■ 10~* 


(3.1 ±0.03) ■ 10"^ 


1.9 ■ lO" ' 


No. of events: Quality > 7 


17 ± 5 


279 ± 3 


204 



TABLE IV: Event numbers for experimental data and Monte Carlo simulations for four major stages in the analysis. The 
errors quoted are statistical only. 



VII. SYSTEMATIC UNCERTAINTIES 

As a novel instrument, AMANDA poses a unique chal- 
lenge of calibration. There are no known natural sources 
of high energy neutrinos, apart from atmospheric neutri- 
nos, whose observation could be used to measure the de- 
tector's response. Understanding the behavior of the de- 
tector is thus a difficult task, dependent partly on labora- 
tory measurements of the individual components, partly 
on observations of artificial light sources embedded in 
the ice, and partly on observations of downgoing muons. 
Even with these measurements, uncertainties in various 
properties that systematically affect the response of the 
detector persist, which prevent us at this time from mak- 
ing a precise measurement of the atmospheric neutrino 
flux. The primary sources of systematic uncertainties, 
and their approximate effects on the number of upgoing 
atmospheric neutrinos in the final data sample, as deter- 
mined by variation of the simulations, are listed below. 

As discusssed in Sections 2 and 3, AMANDA is embed- 
ded in a natural medium, which is the result of millennia 
of climatological history, that has left its mark in the form 
of layers of particulate matter affecting the optical prop- 
erties of the ice. Furthermore, the deployment of optical 
modules requires the melting and refreezing of columns 
of the ice. This cycle results in the formation of bubbles 
in the vicinity of the modules, which increase scattering 
and affect the sensitivity of the optical modules in ways 
that are not yet fully understood. The effects of this 
local hole ice are difficult to separate from the intrinsic 
sensitivity of the OMs. The uncertainties in the neutrino 
rate are approximately 15% from the bulk ice layer mod- 
eling in the Monte Carlo, and as much as 50% from the 
combined effects of the properties of the refrozen hole ice 
close to the OMs, and the intrinsic OM sensitivity, and 
angular response. 

Figure shows two variables that are sensitive to the 
absolute OM sensitivity: the number of OMs hit and the 
velocity of the line fit. The systematic effects of varying 
OM sensitivity on the hit multiplicity for Analysis 1 are 
shown on the top. The peak of the multiplicity distri- 
bution for the standard Monte Carlo (nominal efficiency 
100% — dashed line) lies at a higher value than for the 



data. Reducing the simulated OM sensitivity by 50% re- 
sults in a peak at lower values than the data. The other 
variable strongly affected by the OM sensiti vity - the 
velocity of the line fit, introduced in Section IV A) - is 
the apparent velocity of the observed light front traveling 
through the ice, see Fig. ^ (bottom). 

As a next step, we investigated the effect of the 
angsens OM model (first introduced in Section 3) on 
the atmospheric neutrino Monte Carlo simulation. The 
results of this simulation gave a more consistent descrip- 
tion of the experiment for several variables — e.g., the 
hit multiplicity (the dotted line in Fig. ^) — and they 
produced the absolute neutrino event prediction closer 
to what was found in Analysis 1 (236.9 events predicted, 
223 observed). Similar effects are seen when this Monte 
Carlo is used with Analysis II, however the number of 
predicted events is 25% smaller than observed. Thus the 
cingsens model, while encouraging, does not completely 
predict the properties of observed events in both analy- 
ses. 

Another uncertainty lies in the Monte Carlo routines 
used to propagate muons through the ice and rock sur- 
rounding the detector. A comparison of codes based 
on 1^ and ||ll| indicates that different propagators may 
change the event rates by some 25%. 

Other factors include the simulation of the data ac- 
quisition electronics and possible errors in the time cali- 
brations of individual modules. These effects have been 
studied by systematically varying relevant parameters in 
the Monte Carlo simulations. For realistic levels of vari- 
ation, these effects are well below the 10% level. 

Figure ^ demonstrates how the zenith angle distri- 
bution depends on different atmospheric neutrino event 
generators (our standard generator nusim and an- 
other generator nu2mu |^), and also on the chosen an- 
gular sensitivity of the optical module. Neutrino flavor 
oscillations lead to a further reduction of the angsens 
prediction by 5.4% (in particular, close to the vertical di- 
rection), assuming sin^ 26* = 1 and Am^ = 2.5 • 10"'^ eV^ 
The prediction is reduced by 11% if the largest al- 
lowed Am? is used. 

The combined effect of all these systematic uncertain- 
ties is sufficiently large that simulations of a given at- 
mospheric neutrino flux can produce predictions for the 
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ity level 7 cuts have been applied in all bu t the displayed 



variable (A'^ — 1 cuts, see also Section VP 3 , Fig. 11). The 
vertical dashed lines with the arrow indicate the region of ac- 
ceptance in the displayed variable. In each clear tail 
of fake events is removed by application of the cut, leaving 
good agreement in shape between the remaining events and 
the Monte Carlo expectation. 



0.25 



0.20 



0.15 - 



0.10 



0.05 



0.25 - 



Experiment 
MC: 

OM sensitivity 50% 
OM sensitivity 100% 
angsens OM model 




20 30 40 50 
Number of hit OM (Nch) 



60 



_n Experiment 
MC: 

OM sensitivity 50% 

- OM sensitivity 100% 
angsens OM model 




0.10 0.15 0.20 0.25 0.30 0.35 
Velocity from line fit (vnnefit ) [m/ns] 
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event rate varying by a factor of two. By contrast, the es- 
timated theoretical uncertainty in the atmospheric neu- 
trino flux, at the energies probed by these analyses, is 
30% The effect of neutrino oscillations with the 

Super-K preferred parameters would be less than 10% at 
these energies. 



VIII. SYNTHESIS AND GENERAL OVERVIEW 

Both Analyses I and II are able to separate more than 
200 neutrino event candidates from the 130.1 days of 
AMANDA-BIO detector livetime in 1997. Based on at- 



mospheric neutrino simulations we find that about 4% of 
the total number of events triggered by upward moving 
neutrinos passed the final selection. A total deficit in the 
event rate of about 35% with respect to the standard neu- 
trino Monte Carlo prediction is found for both analyses. 
An event overlap of 102 experimental events is observed, 
consistent with a predicted overlap of 119 ± 13 from the 
atmospheric neutrino Monte Carlo. Thus, the combined 
sample of data provides about 300 neutrino candidates. 
Both analyses estimate their residual background to be 
about 10% of the number of neutrino event candidates. 

Figure ^ shows the energy distribution of the sim- 
ulated neutrinos and the corresponding muon events. 
Ninety percent of all Monte Carlo signal events have 
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FIG. 21: Distribution of simulated zenith angles for different 
neutrino generators and also for a modified angular sensitivity 
of the OM. 



FIG. 23: Effective area for muons versus zenith angle. The 
energy of the muons is given at the point of closest approach 
to the detector. 
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FIG. 22: Energy distributions for simulated atmospheric neu- 
trino events which pass the final neutrino cuts. The effect of 
neutrino oscillations has not been taken into account. The 
figure shows the neutrino and muon energies at the interac- 
tion vertex and the energy of the muons at the point of closest 
approach to the detector center. 



muon (neutrino) energies between 48 (66) GeV and 
1.8 (3.4) TeV. The dominant part of the signal events in 
this analysis comes from neutrino energies below 1 TeV. 
Figure shows the effective area as a function of the 
zenith angle for two ranges of the muon energy at the 
point of closest (POC) approach to the detector. The ef- 
fective area for muons with energies at POC between 100 
and 1000 GeV is 3.9 • 10"^ m^ at trigger level and 2800 m^ 
after application of the neutrino selection cuts. It should 
be noted that much higher effective areas are possible 



when searching for neutrinos from astrophysical point 
sources fs^ or from gamma ray bursts p3| . 

Figure"^ shows the point spread function of the re- 
constructed muon trajectory with respect to the true 
muon direction. Based on Monte Carlo simulations, we 
find a median angular resolution of muons from atmo- 
spheric neutrinos of 3.2° for the final sample. A more 
detailed study of the angular resolution can be found 
in m . Figure ^ shows the skyplot (equatorial 
coordinates) of all the candidate neutrino events found 
across both analyses. The distribution of the events on 
the skyplot is consistent with a random distribution. 

IX. CONCLUSIONS 

The AMANDA-BIO data from 130.1 days of livetime 
during the austral winter of 1997 have been analyzed 
in an effort to detect high energy atmospheric neutrino 
events, and to compare their properties to expectations. 
Two working groups in the collaboration, using differing 
reconstruction, cut optimization and instrumental event 
rejection techniques, produced sets of 223 and 204 neu- 
trino candidates, respectively. Several methods of back- 
ground estimation put the residual event contamination 
from downgoing atmospheric muons and instrumental ar- 
tifacts at about 10%. Taking into account systematic 
uncertainties, the observed event numbers are consistent 
with systematically varied atmospheric neutrino Monte 
Carlo predictions, which are from 150-400 events. The 
range of these predictions is dominated by uncertainties 
in the neutrino flux, in the understanding of photon prop- 
agation through the bulk ice and the refrozen hole ice. 
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and in muon propagation and energy loss. The Monte 
Carlo suggests that 90% of the selected events are pro- 
duced by neutrinos in the energy range of ~ 66 GeV to 
3.4 TeV. The observation of atmospheric neutrinos in 
line with expectations establishes AMANDA-BIO as a 
working neutrino telescope. We finally note that many 
of the procedures for signal separation simplify consider- 
ably in larger detectors. In particular, first results from 
AMANDA-11 1^ demonstrate that the neutrino signal 
is separated with much higher efficiency and with fewer 
cuts than for AMANDA-BIO. 



FIG. 24: Monte Carlo simulation of the angular resolution for 
muons that pass the final selection criteria. The median error 
is 3.2°. 
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FIG. 25: Neutrino skyplot of upgoing events as seen with 
AMANDA-BIO in 1997 in equatorial coordinates. In this fig- 
ure, neutrino events from both analyses are combined. The 
background of non-neutrino events is estimated to be less than 
10%. 
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